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ABSTRACT 

We discuss in detail the broad band X-ray variability expected in Seyfert galaxies in 
the framework of a Compton cooled corona model, computing explicitly spectral indices 
in the 2-10 keV range and fluxes in the medium (2-10 keV), hard (30-100 keV) and 
soft (0.1-0.4, 0.9-2 keV) X-ray bands, for different hypotheses concerning variability 
"modes". In all cases the soft photons responsible for the cooling are assumed to be 
proportionate to the Comptonized flux (fixed geometry), and the plane parallel limit is 
adopted. Variations in the optical depth r of the corona cause spectral changes in the 
Comptonized emission in the sense that the spectrum steepens and the temperature 
decreases with increasing r. If the corona is pair dominated, r is determined by the 
compactness (proportional to luminosity for fixed size). This yields a definite relation 
between spectral shape and intensity implying a variation in intensity of a factor 10 
for a change of 0.2 in the 2-10 keV spectral index. For low pair density coronae, no 
definite prediction is possible. If for instance r varies while the luminosity remains 
approximately constant, the Compton spectrum shows a pivoting behaviour around ~ 
10 keV. In both cases, the coronal temperature decreases for increasing spectral index. 
In the absence of substantial changes in the structure of the corona, the soft thermal 
component correlates with the medium-hard Comptonized component. Although the 
correlation must hold globally, that is for the integrated power of the two components, 
the observed intensity in the soft X-ray band is extremely sensitive to the tempera- 
ture of the thermal component which may vary as a result of variations in luminosity 
and/or area of the reprocessing region(s). Examples are computed using the ROSAT 
response matrix, showing that for kT ;$ 60 eV the spectrum in the ROSAT range soft- 
ens with increasing intensity while the opposite is true for larger temperatures. Current 
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observations of NGC 5548 and Mrk 766 showing different behaviours can therefore be 
reconciled with the model. The spectral behaviour of Mrk 766 strongly suggests that 
the coronal plasma in this source is not pair dominated. The possibility of constraining 
models through simultaneous observations in the soft medium and hard X-ray ranges 
is discussed. 

Subject headings: Accretion Disks - Galaxies: Seyfert - X-Rays: Galaxies 
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1. Introduction 

The observations of OSSE on the Compton 
Gamma Ray Observatory showing a break or an 
exponential cut-off in the high energy spectrum 
of some Seyfert galaxies (Maisack et al. 1993, 
Johnson et al. 1993, Madejski et al. 1995) cou- 
pled with the X-ray observations of GINGA, have 
had a strong impact on our theoretical under- 
standing of the X-ray and 7— ray emission pro- 
cesses in these objects. Though they do not rule 
out pure non-thermal e pair models (Zdziarski, 
Lightman & Maciolek-Niedzwiecki 1993), ther- 
mal or quasi-thermal Comptonization models 
are favoured (e.g. Sunyaev & Titarchuck 1980; 
Haardt & Maraschi 1991 and 1993, hereinafter 
HM91 and HM93; Ghisellini, Haardt & Fabian 
1993; Titarchuck & Mastichiadis 1994; Zdziarski 
et al. 1994; Zdziarski et al. 1995; see Svensson 
1996 for a recent review). 

The origin of the Comptonizing electrons is 
at present unknown. They could be related to 
an accretion disk corona heated by magnetic 
dissipation processes and/or they may be elec- 
tron positron pairs produced by photon photon 
collisions (see e.g. Fabian 1994 and references 
therein). The electron temperature can be con- 
strained by a thermal balance equation, where 
the heating rate is the total (direct plus repro- 
cessed) luminosity and the cooling mechanism is 
Comptonization of soft photons. 

The hot corona is probably coupled to a cooler 
optically thick gas layer (presumably the accre- 
tion disk itself), which has three major roles: i) it 
reprocesses a fixed fraction of the Comptonized 
radiation into a soft blackbody component; ii) a 
fixed fraction of the reprocessed blackbody pho- 
tons crosses the hot region acting as seed pho- 
tons for the Comptonization process; iii) it re- 
flects medium energy X-rays giving rise to the 
so called reflection component and Fe line emis- 
sion. This picture (HM91, HM93) is suggested 
by the presence of the Compton reflection hump 
(Pounds et al. 1990) and by the relatively small 



dispersion in the distribution of spectral indices 
of Seyfert I galaxies (Nandra & Pounds 1994). In 
fact, if X-ray reprocessing contributes soft pho- 
tons proportionately to the medium-hard X-ray 
luminosity, the 'feed-back' effect keeps a quasi- 
constant X-ray spectral shape, weakly dependent 
on the source luminosity. The broad band spec- 
tra derived from this model are consistent with 
the average Seyfert spectrum derived from OSSE 
observations (Zdziarski et al. 1995; Gondek et al. 
1996). 

It is fair to point out that the best studied 
object, NGC 4151, shows (at least in low state) 
an unusually flat X-ray spectral index (a ~ 0.5) 
and a low temperature (kT ~ 60 keV) imply- 
ing a Comptonized-to-thermal luminosity ratio 
L c /L s ^ 10, far from the limit L c /L s ~ 2 ob- 
tained in a plane parallel geometry. This may 
be accounted for by a special geometry, for in- 
stance if the high energy emission originates in 
localized non planar active blobs as in the model 
proposed by Haardt, Maraschi & Ghisellini 1994 
(hereinafter HMG). Alternative models involv- 
ing a mirror reflection of the observed radiation 
(Poutanen et al. 1996) or reprocessing in clouds 
along the line of sight (Zdziarski & Magdziarz 
1996) may also apply for this special object. 

In Comptonization models the spectral shape 
in the X-ray range is mainly determined by the 
combination of the temperature O = kT/mc 2 , 
and the optical depth of the scattering electrons 
r, while the cut off energy is related essentially 
to O. Simultaneous measurements of the spectral 
index a and of the cut off energy can therefore 
determine the physical parameters of the Comp- 
tonizing region^. 

"Spectral dynamics", i.e., the variation of the 

lr The exact relation between a and r depends on the ge- 
ometry of the scattering region. The optical depth can 
therefore be deducted from observations only "modulo ge- 
ometry". On the other hand fitting the data assuming a 
simple homogeneous sphere gives a value of r that is also 
a good estimate of the average (over all directions) scat- 
tering optical depth for any other possible geometry. 
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spectral shape with intensity is an additional ba- 
sic diagnostic tool (see, e.g., Netzer, Turner &; 
George 1994) which can test the model and pro- 
vide insight into the way the corona is heated 
and cooled. Some of the questions that could be 
answered are: is the optical depth of the corona 
dominated by e + e~ pairs? How do the optical 
depth and electron temperature vary with lu- 
minosity? Are the reprocessed components pro- 
portionate to the Comptonized one during vari- 
ations? For instance a careful analysis of EX- 
OSAT spectra of Seyfert galaxies, showed that 
the observed spectral variability is difficult to 
reconcile with non-thermal pair models (Grandi, 
Done & Urry, 1994) . 

Another important prediction of the model 
is the presence of a thermal component which 
should be proportionate to the Comptonized one 
in the limit of a constant geometry. Its temper- 
ature however is not completely constrained by 
the model and could range from UV to soft X- 
rays. The presence of a "soft excess" in the X-ray 
spectra of AGN has been detected by different ex- 
periments (EXOSAT, ROSAT and more recently 
ASCA) (e.g. Arnaud et al. 1985; Pounds et 
al. 1986; Wilkes k Elvis 1987; Turner & Pounds 
1989; Elvis, Wilkes & Mc Dowell 1991; Comas- 
tri et al. 1991; Saxton et al. 1993; Walter & 
Fink 1993; Turner, George & Mushotzky 1993; 
Brandt et al. 1994; Done et. al 1995). It has 
been interpreted in some cases as a genuine ther- 
mal component possibly due to reprocessing (e.g. 
Leighly et al. 1996), or as the high energy tail of 
the UV bump (e.g. Walter & Fink 1993), while 
in others as an artefact of the residuals due to 
incomplete absorption (warm absorber, see e.g. 
Done et al. 1995). Instruments with improved 
energy resolution in the soft band are needed 
to clarify the interpretation of the "soft excess" 
which may even involve more than one compo- 
nent. However a distinctive point between differ- 
ent models is the predicted correlation between 
the soft and medium X-ray intensities. On this 
aspect useful data obtained from extensive ob- 



servations of some Seyfert galaxies with ROSAT 
and/or ASCA already exist and more will surely 
become available in the near future. 

We are therefore motivated to examine in de- 
tail the broad band spectral variability expected 
in the accretion disk with hot corona model de- 
veloped in previous papers (HM91, HM93 and 
HMG). We will concentrate on the case of extra- 
galactic objects such as Seyfert I galaxies. The is- 
sue is complicated for several reasons. A quanti- 
tative analysis requires spectral simulations since 
"second order effects" are important (see §3) . 
However the main uncertainty is our poor knowl- 
edge of the physics of the corona. For instance 
do the active coronal regions vary in size and/or 
height above the reprocessing disk? And how is 
the heating rate related to the optical depth? 

In the following we chose the simple case in 
which the active and reprocessing region(s) can 
vary in area but maintain a fixed ratio between 
Comptonized and soft cooling photons. In par- 
ticular we choose the limit of a plane-parallel 
geometry characterized by a Comptonized-to- 
soft ratio L c /L s ~ 2. This is not inconsistent 
with an inhomogeneous picture where each ac- 
tive region can be approximated as plane paral- 
lel. With these restrictions we compute explic- 
itly the expected relations between "observable 
quantities", e.g. the variation of the 2-10 keV 
spectral index (including reflection) with temper- 
ature and intensity, the variation of the soft in- 
tensity "S" (0.1-2 keV) with medium "M" (2-10 
keV) and hard "H" (E > 30 keV) intensity. 

The plan of the paper is as follows. In §2 
we describe the model and the methods used 
for the computation of the spectra and for the 
derivation of the "observable" quantities. Since 
the key parameter driving spectral variability in 
the medium band is r, while spectral variabil- 
ity in the soft band is essentially related to Tbb, 
we found convenient to discuss relations between 
spectral and flux variability in the two bands sep- 
arately. For this reason, in §3 we deal with the 
expected spectral variability in the medium and 
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hard X-ray bands, using the 2-10 keV spectral 
index, the 2-10 keV intensity, the temperature 
of the hot corona and the 30-100 keV intensity, 
while in §4 we compute the intensity and spec- 
tral variations of the soft reprocessed blackbody 
component for different variations of the system 
parameters. The results are converted into count 
rates and hardness ratios in the 0.1-0.4 and 0.9- 
2 keV bands using the ROSAT response matrix 
and a standard neutral hydrogen column. We 
discuss the expected correlations between vari- 
ability in the soft and medium energy bands and 
compare the results to available observations in 
§5. An exhaustive summary and final conclusions 
are presented in §6. 

2. The Disk Corona Model 

We follow a scenario originally proposed by 
Liang (1979), and revived by HM91 and HM93, 
in which the accretion flow is illuminated by 
the hard X-ray source (as the observed reflec- 
tion spectrum indicates). A large fraction (about 
80—90 per cent) of the incident power is expected 
to be thermalized in the disk and re-emitted as 
thermal reprocessed radiation, which is in turn 
Comptonized in the corona producing the X- 
rays. If the reprocessed radiation energy density 
dominates over the local disk thermal emission, 
the Compton y parameter is forced to be ~ 1 
independently of the coronal scattering opacity. 
This in turn ensures that the Comptonized radi- 
ation exhibits a spectral index close to 1. One 
possible configuration which satisfies the above 
condition is a plane parallel corona above an ac- 
cretion disk. In this case the condition y ~ 1 is 
achieved only if the entire available gravitational 
power is released in the hot corona. In this case 
the reprocessed thermalized radiation coincides 
with the UV bump, and the reprocessed reflected 
radiation forms the 30 keV Compton hump. 

The same condition can be satisfied if there 
are several such smaller regions above the disk, 
rather than a single smooth corona (HMG). In 
the latter case only a fraction of the available 



energy needs to be released in the hot corona. 
However being localized in hot spots, the repro- 
cessed thermalized radiation can still dominate 
the cooling. The underlying accretion disk can 
contribute to most or at least part of the total 
luminosity in the UV band as in the standard 
Shakura-Sunyaev disk model (Shakura &; Sun- 
yaev 1973). As in the uniform corona case, the 
reprocessed reflected radiation forms the 30-keV 
Compton hump, while the reprocessed thermal- 
ized radiation is probably emitted in hot spots in 
the EUV-soft X-ray band (HMG). Therefore the 
emission in the EUV-soft X-ray band should be 
closely connected to the hard X-ray radiation. 

In the present paper we assume that the cool- 
ing time of the electrons in the corona is shorter 
than the time scale of variation of the heating 
rate, so that we can consider spectral and flux 
variability as a succession of stationary states. 
The Comptonized spectra are computed in plane 
parallel geometry for different inclination angles 
A quadrature scheme was implemented to 
solve the radiation field for successive scattering 
orders, taking into account anisotropic Comp- 
ton emission for the first scattering distribution 
(Haardt 1993). The single scattering photon dis- 
tribution is computed by means of the fully rel- 
ativistic kernel for isotropic unpolarized radia- 
tion (Coppi & Blandford 1990). This ensures 
that the exponential roll over of the spectrum 
is correctly modeled as we tested against Monte- 
carlo calculations. The present numerical spectra 
are therefore in much closer agreement with re- 
cent detailed calculations by Poutanen & Svens- 
son (1996) than the approximated treatment of 
HM93. An extensive discussion of the method 
we used can be found in Haardt (1994). 

The downward Compton flux (i.e. the X-ray 
flux leaving the corona with an angle •& > vr/2) is 
averaged over the solid angle, and then reflected 
by an infinite neutral layer (which mimics the ac- 
cretion disk) with solar metal abundances. The 
angular distribution of the reflected spectrum is 
computed in the approximated way described in 
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Ghisellini, Haardt & Matt (1994) (see Burigana 
1995 and Magdziarz & Zdziarski 1995 for detailed 
calculations). The angular dependent reflected 
spectrum is then Comptonized in the corona and 
added to the primary continuum. The "measur- 
able" spectral index in the [2-10] keV band is 
obtained fitting with a power law to the com- 
posite spectrum disregarding the goodness of the 
fit. The photon flux (i.e. number of photons per 
unit time per unit surface) in interesting bands 
is readily obtained by integration over energy. 

Concerning the soft photon input, we con- 
sider cases in which the reprocessed radiation is 
isothermal and has a temperature Tb b comprised 
between 10 and 300 eV. This scenario can be re- 
garded as a uniform corona located above the 
very inner part of the accretion disk, or alter- 
natively as a blobby structure where the overall 
X-ray emission is dominated by few blobs with 
similar properties. 

With our approximations the interesting quan- 
tities are computed with good accuracy in a quick 
way. A more detailed treatment of the full radia- 
tive transfer problem in the same context, includ- 
ing free-free radiation, double Compton scatter- 
ing, detailed pair balance and polarized scat- 
tering, can be found in Poutanen & Svensson 
(1996). 

3. Variability in the Medium and Hard 
X Ray Domains 

3.1. The a — t Relation 

Since the model constrains the Compton pa- 
rameter to be almost constant, one expects that 
the spectral index of the Comptonized spectrum 
is fairly constant too, and does not depend on 
the three parameters of the system (i.e. r, 6 
and T B b)- 

The above argument is however valid only as 
a first order approximation, in particular in the 
case of a disk-corona system, where the scatter- 
ing anisotropy is important and when we consider 
the actually observed spectrum which includes 



the reflection hump. We therefore computed the 
spectral index of the total spectrum in the [2-10] 
keV range ccp-ro] a s described in the previous 
section for different values of the coronal param- 
eters and of the inclination angle. Comptoniza- 
tion of the reflected spectrum is also included. 

The results are shown in Fig. 1 and Fig. 2 
where ap-io] is plotted vs. the optical depth 
and the temperature of the corona respectively, 
for two values of the inclination angle and two 
values of the temperature of the soft photons. In 
order to show the effect of reflection on the total 
spectrum, in Fig. 3 we plotted, as a function of r, 
the difference between the spectral index of the 
intrinsic continuum alone and that of the total 
spectrum (continuum+reflection) . 

The optical depth was varied between 0.1 and 
1, corresponding to an equilibrium temperature 
of the hot corona in the range 30-300 keV, as 
suggested by observations (Maisak et al. 1993, 
Johnson et al. 1993, Madejski et al. 1995). 

The important points to note are the follow- 
ing: 

1. Changes in r give rise to significant spec- 
tral variability despite the fact that the ra- 
tio L c /L s is constant in the model. In the 
plane parallel limit considered here ap-io] 
is found to be in the range [0.4-1.4] for r 
varying between 0.1 and 1. 

2. Steeper spectra are produced as r increases 
approaching unity. 

3. An increase of the spectral index is accom- 
panied by a decrease of the coronal tem- 
perature, i.e. spectral index and coronal 
temperature are anticorrelated. The ac- 
tual values of a and kT in Fig 2 refer to 
L c /L s ~ 2, however the anticorrelation 
between spectral index and the electron 
temperature is a general feature common 
to models based on Compton cooling with 
fixed L c /L s . For higher values of L c /L s a 
similar curve would be obtained, but with 
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values of a shifted downwards. As shown 
in Fig. 2, the dispersion of the ap-io] vs - 
kT relation is larger at higher temperature. 

A physical understanding of these results in- 
volves several effects: 

1. for any given r the temperature has to 
adjust to maintain an almost fixed ratio 
L c /L s , and for increasing r the tempera- 
ture decreases. The resulting average spec- 
tral index (i.e. the spectral index averaged 
over the total solid angle) is, for the plane 
parallel geometry adopted here, slightly 
larger than 1 for any value of r (HM93). 
Comptonization theory implies that, in or- 
der to keep a constant ratio L c /L s , a spec- 
tral index steeper than 1 must increase with 
decreasing (i.e. increasing r). This can 
be understood noticing that the extrapola- 
tion to low energy of the Compton power 
law must intercept the peak of the injected 
black body soft photon distribution. Since 
for a > 1 the integral over energy of the 
power law spectrum mainly depends on the 
lower integration limit E\ oc (160 2 + 40 + 
l)kTBB, a shift to lower energy of E\ (i.e. a 
decrease of 0) keeping a constant a would 
result in a larger ratio L c /L s . Therefore, in 
order to keep the ratio constant, the spec- 
tral index is forced to steepen (for the same 
reason, if a < 1, a reduction of implies 
a further flattening, rather than a steep- 
ening, of a). The steepening of a halts 
when a further decrease of does not cause 
a shift of Ei to lower energy, i.e. when 
160 2 + 40 < 1, say 160 2 + 40 < 0.5. For 
even smaller temperatures the spectral in- 
dex must flatten (and hence there is a max- 
imum in the a vs. curve) since, though 
Ei is constant, the upper integration limit 
(oc 0) decreases anyway. In terms of r, 
since in our adopted geometry the Comp- 
ton parameter y ~ (160 2 + 40)r is ~ 0.6 
(HM93, Poutanen & Svensson 1996), the 



"turning point" occurs for r ~ 1. Note that 
for a < 1 {L c /L s ^ 2) there is a monotonic 
flattening of a as decreases. The effect 
discussed here is independent of the par- 
ticular geometry adopted. In brief, a fixed 
L c /L s ratio does not imply an exactly con- 
stant spectral index [e.g. see Fig.2 in Ghis- 
ellini & Haardt (1994; hereinafter GH94) 
for the simplest possible geometry, i.e. a 
homogeneous spherical source, where this 
effect is present]; 

2. the outgoing Compton spectrum in plane 
parallel geometry is a function of the emis- 
sion angle, because of the Compton rocket 
effect (Haardt 1993). Flatter spectra are 
seen for lower inclinations, and the impor- 
tance of the anisotropic Compton emission 
decreases with decreasing coronal tempera- 
tures. Anisotropic Compton emission is the 
main cause of the reduction of ap-io] f° r 
small r. This effect goes in the same direc- 
tion as that discussed in the above point 1, 
but, contrary to it, depends on the adopted 
geometry (e.g. it is absent in GH94). The 
larger dispersion of the ccp-io] vs - kT rela- 
tion at high temperature (see Fig. 2) is due 
to the increasing importance of anisotropic 
Compton emission for increasing kT (see 
Haardt 1993 for details); 

3. the contribution of the reflected Compton 
hump to the emitted spectrum is oc e~ T ^, 
where \i is the cosine of the inclination 
angle, so it is smaller for larger r (for 
fixed inclination angle), and can be boosted 
by anisotropic Compton scattering (HM93) 
for high temperatures. As long as r £ 0.4, 
the contribution of reflection to the [2-10] 
spectrum is marginal (Aa^-io] — 0-05 — 
0.08). For lower r the electron tempera- 
ture is high enough for anisotropic Comp- 
tonization to boost reflection. The effect, 
for fixed r, is larger for larger fi and larger 
kTsB (see Fig. 3). Again, this is an ef- 
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feet whose net result is to flatten the total 
spectrum as r decreases; 

4. finally we discuss the effect on a of a change 
in Tbb- In fact the Compton y param- 
eter depends not only on L c /L s but also 
(only weakly when y 1) on the mean en- 
ergy of the soft photons. This can be eas- 
ily understood considering that the higher 
the energy of the soft photons injected in 
the corona, the lower the number of scat- 
terings needed to reach an energy ~ kT. 
When y <C 1 the (low) cooling is due only 
to the first scattering, so that L c /L s ~ y. 
The relation between L c /L s and y does not 
depend on any other parameter. But as 
soon as y increases higher order scatter- 
ings become important in the cooling pro- 
cess. The asymptotic limit for y ^> 1, 
i.e. when all the photons reached ther- 
mal equilibrium, is L c /L s = 40 /xq where 
xq is the dimensionless soft photon mean 
energy ~ 2.7kTBB /mc 2 . For this reason, 
for a given r the equilibrium temperature 
necessary to sustain a given L c /L s is (only 
slightly for y ~ 1) higher (and hence the 
Compton spectrum flatter) for higher Tbb 
(for a visualization of how y varies with 
Tbb see Fig. 1 in Shapiro, Lightman k, 
Eardly 1976). Large chang es in Tb b j un- 
likely to occur in a given source, are re- 
quired to produce significant changes in a. 
However, this effect could be relevant when 
comparing different classes of sources (e.g. 
Seyferts and QSOs) where the difference 
between the average values of Tbb can be 
large enough to influence the value of the 
spectral index of the hard X-rays. As we 
will discuss in §4, variations of Tbb m a 
particular object, even if small, can play 
an important role in the observed variabil- 
ity in the soft X-ray band. 

As discussed in §2, ap-io] is a simple least 
square fit to the total spectrum, which is not a 



perfect power law (because of, e.g., anisotropy 
and inclusion of reflection). In order to quan- 
tify the differences between the computed spec- 
tra and a power law, we associated a fixed per- 
centual error to each of our 10 computed spectral 
points. By varying the error until a reduced x 2 
of unity was obtained in the least square fits, we 
estimated a minimum observational precision re- 
quired to detect the spectral distortions. For the 
most distorted spectra (those at low inclination 
and low r) an observational error £,2% would 
result in a reduced x 2 larger than unity. The 
typical maximum error allowed is ^1%. 

We now discuss the relation between ap-10] 
and the X-ray luminosity L c on the basis of two 
alternative limits. The first case corresponds to 
a pair dominated corona where the optical depth 
is determined by the compactness. In the second 
case the relation between r and L c is essentially 
unknown, and we will consider the limit case in 
which r varies with constant luminosity. 

3.2. The a— Flux Relation in a Pair Dom- 
inated Corona 

If the corona is pair dominated, as plausi- 
ble in the case of compactnesses i c ^ 10, where 
l c = (L c / R)(<jt /m e c 3 ) and R is the typical size- 
scale of the emitting region (see Svensson 1996 
for an exact definition of compactness in different 
geometries), the optical depth of the hot phase 
are determined by the compactness alone. As 
discussed in HM91, once the ratio L c /L s and 
the absolute value of i c are specified, there exists 
a univocal correspondence of these parameters 
with r, under the assumption that pairs domi- 
nate the scattering opacity over normal plasma. 
Thus if variations of the compactness are due to 
variations of the dissipated luminosity L c , the op- 
tical depth and hence the actual value of a [2- 10] > 
depend on the intensity of the source]^] 

2 Stern et al. (1994) and more recently Poutanen & Svens- 
son (1996) and Svensson (1996) pointed out that the ex- 
act relation between the rate of pair production and the 
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From the ap-io] VS- r elation computed pre- 
viously (Fig. 1) we can therefore derive for a 
pair dominated corona an ap-io] vs - intensity 
relation. This is shown in Fig. 4 for two val- 
ues of the viewing angle. It is quite clear (see 
also Fig. 2 in Stern et al. 1996) that noticeable 
variations of ap-io] ( sa Y ^ a [2-io] — 0-3) occur 
as the luminosity (and hence the compactness) 
varies by at least a factor of 20. Thus for a pair 
dominated corona the spectral index must be es- 
sentially constant during "normal" AGN X-ray 
variability (i.e. flux variations within a factor of 
2 or so). Modest spectral steepening can occur 
during very large intensity increases. 

On the other hand a spectral steepening also 
implies a decrease of the coronal temperature 
(Fig. 2) (see also Pietrini h Krolik 1995). Since 
an observational determination of the latter with 
presently available instrumentation is not easy, 
we also computed the variation in the hard pho- 
ton intensity ([30-100] keV), that would be as- 
sociated with a variation in the [2-10] keV in- 
tensity. This is shown in Fig. 5, where we also 
report values of ap-io] corresponding to different 
values of the intensity. The two intensity scales 
are arbitrary but their ratio is not. Due to the 
temperature decrease for increasing luminosity, 
the intensity in the hard band varies less than 
that in the medium band. However the two are 
positively correlated. Because of the large over- 
all luminosity variation the steeper and cooler 
spectrum does not cross the flatter, hotter one, 
except at very high energy. 

Different conclusions hold if variations of i c 
are due to variations of the linear dimension of 
the source R rather than to variations of the in- 
trinsic luminosity. Such a circumstance can oc- 
cur if the corona is formed by many active regions 
rather than being homogeneous. In this case sig- 
nificant spectral variations could occur without 

heating/cooling rates depends on the details of the calcu- 
lations. We have taken the compactness vs. r curve given 
by Svensson (1996) in order to compute the photon flux 
in the [2-10] keV band. 



large variations in the output. Strong effects on 
the soft X-ray spectrum are expected since the 
temperature of the reprocessed radiation scales 
as 1/R 2 . A detailed discussion will be presented 
in the next §4. 

3.3. The a— Flux Relation in a Low Pair 
Density Corona 

If the compactness is ^ 10, the corresponding 
pair optical depth is ;$ 0.1 (Stern et al. 1994). 
The main contribution to the optical depth then 
comes from "normal" electrons, whose amount is 
not related to the source luminosity in a simple 
and easily predictable way. We therefore consider 
the case in which the optical depth changes while 
the dissipated luminosity does not vary substan- 
tially. 

In Fig. 6 we show the total spectra for some 
choices of input parameters. It is evident from 
the figure that, despite the fact that the total 
luminosity in the spectra is the same (when in- 
tegrated over all the emission angles) , the flux in 
the [2-10] keV band varies. The spectral index 
vs. the photon flux in the [2-10] keV band is 
shown in Fig. 7. We can see that small observed 
variations of the count rate can be accompanied 
by significant spectral variations. Another im- 
portant result, only weakly dependent on details, 
is that as long as ap-io] ~ 1, it correlates to the 
photon flux, i.e. the spectrum softens with in- 
creasing intensity, while for steeper spectra the 
two quantities are anticorrelated. 

From Fig. 6 it is also evident that the flux 
above 50 keV is going to decrease for increas- 
ing spectral index. The decrement is relatively 
modest as long as ap-io] ~ 1 since most of the 
power is carried by the few high energy photons, 
but becomes large as ap-io] ~ 1- Combining the 
spectral index intensity relation shown in Fig. 
7 with the spectral index temperature relation 
we can compute the intensity in the hard X-ray 
band ([30-100] keV) vs. the [2-10] keV inten- 
sity. This is shown in Fig. 8. For ap-io] ~ 1 
the emission in the [2-10] keV band is positively 
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correlated with the emission in the 7— ray band: 
small increments in the X-ray flux are accom- 
panied with larger increments in the 7— ray flux. 
For ap-io] ~ 1 the correlation is in the opposite 
sense, and is less strong: small increments in the 
medium X-ray emission relate to small decre- 
ments of the flux in the 7— ray band. Around 
the turning point a ~ 1 we have the largest vari- 
ation of the 7-ray flux for the smallest variation 
of the X-ray flux. 

In brief, the model predicts that if noticeable 
spectral variability is observed while the luminos- 
ity variation is less than factor of ~ 2, then the 
scattering optical depth of the corona is not ( or 
weakly) related to the intrinsic coronal luminos- 
ity. This may lead to two possible conclusions: 

1. pairs are not important as a source of scat- 
tering material, 

or 

2. if pairs dominate the opacity, the linear 
dimensions of the source must vary, i.e. 
the compactness varies while the luminos- 
ity stays approximately constant. 

In the latter case the temperature of the soft 
reprocessed photons must vary, increasing when 
the medium energy X-ray spectrum steepens. In 
any case, strong spectral variability can well oc- 
cur during "normal" (within a factor of 2) flux 
variations. 

We caution here that in order to measure a 
variation in luminosity one should measure the 
intensity up to the high frequency cut off, that is 
up to the 100 keV range. In the absence of such 
information one can estimate the luminosity vari- 
ation from the known medium energy intensity 
and the expected high energy cut off. 

4. Variability in the Soft X Ray Domain 
and Correlations to the Medium Band 

One of the basic features of the model is that 
a fraction of the Comptonized X-rays is repro- 



cessed into soft thermal photons which are in 
turn responsible for the cooling of the corona. 
This fraction is fixed for fixed geometry, therefore 
the model predicts a "perfect" correlation be- 
tween the Comptonized luminosity and the soft 
reprocessed one. Observationally, the latter com- 
ponent could be identified with the soft excess 
observed in Seyfert galaxies, although there is 
no general consensus on the very nature of the 
soft component. The soft X-ray emission has 
been identified as the high energy tail of the 
UV bump (Walter & Fink 1993), but such claim 
and in general whether or not the soft X-rays 
are correlated to the UV emission is still contro- 
versial (see, e.g., Laor et al. 1994 and Ulrich- 
Demoulin & Molendi 1996a, 1996b). The low 
ROSAT PSPC energy resolution does not allow 
to constrain the spectral shape of the excess, and 
usually a steep power law is fairly adequate to fit 
the data (Gondhalekar, Rouillon-Foley & Kel- 
let 1996). Higher resolution ASCA data seem to 
support a thermal origin of this component (e.g. 
Leighly et al. 1996), though absorption features 
due to highly ionized gas along the AGN line of 
sight could be responsible (at least of part) of 
the excess of residuals in the single power law 
fits (e.g. Cappi et al. 1996). Indeed some re- 
cent ASCA data are consistent either with an 
emission thermal component or with absorption 
due to highly ionized gas (IC 4329A, Cappi et 
al. 1996; NGC 3227, Ptak et. al. 1995). In sev- 
eral objects both a warm absorber and a thermal 
emission component seem to be present (NGC 
5548, Done et al. 1995; NGC 4051, Guainazzi et 
al. 1995; see also Netzer, Turner k, George 1994). 
Mathur, Elvis & Wilkes (1995) argued that a 
highly ionized, high velocity outflowing gas re- 
sponsible for X-ray and UV absorption might be 
a common component in Quasars and should be 
located outside the Broad Emission Line Region. 

At least for NGC 5548 there are indications 
that the soft emission and the power law vary in a 
correlated fashion, but the correlation is non lin- 
ear, i.e. the soft component varies more (Nandra 
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et al. 1993, Done et al. 1995). In one case a 30% 
increase of the soft counts observed with ROSAT 
was not accompanied by a significant variation in 
the estimated power law component. This event 
suggested that at least soft X-ray "flares" may 
well have a different origin than reprocessing of 
hard x-rays (Done et al. 1995). 

We wish to stress here that the behaviour 
of the observed soft X-ray intensity depends 
strongly on the spectrum, i.e. the temperature 
of the reprocessed emission Tbb- The latter de- 
pends on the absolute luminosity of the X-ray 
emission (we recall that in the model it equals 
roughly half of the total power produced via 
Comptonization) and on the size of the radiating 
surface. Unfortunately we do not know whether 
and how these two quantities are related with 
each other. As a first guess we may expect that 
both R and L scale with the mass of the source. 
This may determine a mean value of L/R? for a 
given source ( sources with smaller masses would 
have hotter soft components), however how the 
ratio varies when the luminosity varies in a given 
source remains highly unpredictable. Further- 
more, the amplitude of variations in fixed energy 
bands can be different e.g. the hardness ratios 
can change substantially when the response of 
a particular detector is considered. Comparison 
with observations can then be made only "filter- 
ing" simulated light curves through the actual 
response matrix of existing experiments. 

In brief, although there is a one-to-one cor- 
respondence between the Comptonized and re- 
processed emission, this holds on global scale, 
i.e. considering the integrated total luminosities. 
What is actually observed as "count rate" in the 
given bands can be very different from what is ex- 
pected on this simple basis. In the following we 
compute the expected behaviour of the soft X-ray 
intensity in response to variations of the medium 
energy intensity in two limiting cases, for vary- 
ing luminosity and fixed size of the reprocess- 
ing region and for constant luminosity and vary- 
ing size. Since most of the observational data 



presently available derive from ROSAT we ex- 
plicitly compute intensities and hardness ratios 
in the most commonly used ROSAT bands. 

4.1. Simulated ROSAT light curves 

We performed simulations of variability pat- 
terns according to selected rules that will be de- 
scribed in the following. The time dependent 
broad band spectra were then filtered through 
a Galactic neutral absorber assuming Nh = 
2 x 10 20 cm" 2 , and finally folded with the ROSAT 
response matrix. We divided the ROSAT band 
in three intervals, defining a "hard count rate" as 
the count rate in the [0.9-2] keV band (C[ .g_2]; 
it should be kept in mind that in all the other 
sections the "hard" band refers to energies ^ 30 
keV) and a "soft count rate" as the count rate 
in the [0.1-0.4] keV band (C[ .i-o.4])- We then 
computed the hardness ratio defined as HR= 
C[o.9-2]/C[o.i-o.4] • This quantity provides direct 
spectral information in the ROSAT range. 

It is important to realize that the dependence 
of HR on the temperature of the soft component 
is not monotonic. This is shown in Fig. 9. a for 
two cases with regard to the power law compo- 
nent (both with 0° inclination), one correspond- 
ing to low optical depth r = 0.1 (resulting in 
a flat spectral index a^-io] —0.5-0.6, see Fig. 
1), the second to higher optical depth r = 0.65 
(yielding to ap-io] ~1. 2-1.3) For a pair domi- 
nated source these values of r would imply a large 
difference in compactness, roughly 10 3 larger for 
the high r case (see Svensson 1996). Interme- 
diate cases must be comprised between the two 
curves. The relation is independent of the actual 
mechanism causing the variations of Tg b , but de- 
pends somewhat on the parameters defining the 
hard power law. 

The behaviour of HR vs. UTbb is easily un- 
derstandable. For low temperatures, only the 
exponential tail of the Planck distribution falls 
in the ROSAT band. In this regime a temper- 
ature increase leads to a large increase of the 
count rate in the soft band only: the spectrum 
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becomes softer. The softening stops when the 
peak of the blackbody distribution is within the 
band, and the exponential tail starts to affect the 
intensity in the [0.9-2] keV band. In the latter 
regime when the temperature increases the spec- 
trum becomes harder. The transition occurs for 
kTsB around 60 eV. 

In order to compute the relation between Tbb 
and the count rate in the ROSAT band, we have 
to specify the mechanisms leading to the tem- 
perature variations. We have considered two dif- 
ferent ideal ways in which Tbb can vary, the 
first corresponding to variations in luminosity for 
fixed size, the second to variations in size for fixed 
luminosity, and discuss them in turn. 

In all the cases discussed in the present sec- 
tion the reprocessed thermal component is su- 
perimposed to a standard multicolor disk emis- 
sion arising from a 1O 7 M black hole radiating 
at 10% of the Eddington limit. The X-ray lu- 
minosity is taken as 0.3 of the disk emission and 
a face-on line of sight is assumed. We checked 
that our results are only marginally sensitive to 
disk parameters. For rather extreme cases such 
as a 1O 6 M black hole accretion disk radiating at 
the Eddington limit, with an X-ray luminosity of 
only 10% of the UV disk emission, the ROSAT 
band is still dominated by the reprocessed X-rays 
as long as the temperature of the latter compo- 
nent is £ 20 eV. The net effect of increasing the 
inclination angle is to reduce the amplitude of 
variations of the count rate in the soft band. 

4.2. Variations of Intrinsic Luminosity 

Suppose the luminosity of the source in the 
power law component changes (i.e. L c changes), 
while r and remain constant. In this case 
an increase in luminosity produces an increase 
in temperature of the thermal component, and 
hence a correlation between medium and soft X- 
rays is expected. 

The hardness ratio plotted in Fig. 9. a can now 
be plotted as a function of the intensity which 



would give rise to the considered temperature 
variation. This requires a very large intensity 
range, a factor 10 for kTsB in the range 10-316 
eV. , implausible for the same source. Thus a 
single source will be confined to some restricted 
interval in intensity, while different sources could 
fall in different regions of the plot. The results 
are shown in Fig. 9.b. The two curves correspond 
to the same coronal parameters used in panel (a). 
Clearly the count rate scale is arbitrary, but the 
absolute value of the hardness ratio is meaning- 
ful, being independent of the source luminosity 
for a given spectral shape. For reference, we plot 
along each curve discrete points corresponding to 
an increase in luminosity by a factor 2, i.e. to an 
increase of 20% in Tbb with respect to the pre- 
ceding one. The exact value of VTbb relative to 
each point can be easily read from Fig. 9. a where 
the same discrete points are marked. 

The count rates in the two bands are always 
correlated, though the correlation is not linear. 
Even though the temperature only varies as 
as long as the temperature is ^ 60 eV the num- 
ber of soft photons emitted in the ROSAT band 
changes significantly for even modest variations 
of temperature. This regime corresponds to the 
decreasing branches of the curves in Fig. 9.b. 
Along this branch an increase in luminosity pro- 
duces an exponential increase of C[ .i_o.4], while 
the increase of C[ .g_2] is linear, being dominated 
by power law photons. The hardness ratio is 
clearly anticorrelated to the count rate. The to- 
tal count rate shows large variations, being dom- 
inated by the lowest energy channels. 

The behaviour is very different when the peak 
of the black body spectrum is well within the 
ROSAT band, i.e. for kT > 60 eV. The hardness 
ratio is constant as long as the exponential tail of 
the Planck distribution does not affect the [0.9-2] 
band. When this happens, C[o.9-2] responds ex- 
ponentially to AL, while C[ .i_o.4] varies almost 
linearly. The total count rate is almost linear 
with AL, so the spectrum becomes harder. 

The case of a steep power law component is 
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qualitatively similar but less extreme. The vari- 
ations of the hardness ratio are smaller in both 
branches. 

In a pair dominated source luminosity varia- 
tions with constant coronal parameters are, strictly 
speaking, not allowed. States with different lumi- 
nosity belong to curves with different values of r. 
However, since the separation of the two curves 
for constant Tbb (dotted lines) corresponds to 
an increase of £ c of a factor ~ 10 3 , the shift cor- 
responding to a factor 2 is clearly negligible. 

4.3. Variations of Linear Dimensions 

We consider the case of a constant output ra- 
diating from a varying surface area, while r and 
remain constant. This may represent a case 
in which the primary hard emission occurs at a 
constant rate in an expanding (and/or contract- 
ing) corona, or blob, or if the emission occurs 
at different time in blobs with roughly the same 
luminosity but different sizes. 

Though the luminosity of the hard and soft 
components is constant, the temperature of the 
thermal photons will change giving rise to observ- 
able effects in the soft X-ray domain. To study 
the problem we varied the linear size R of the 
emitting region, so that the resulting black body 
temperature covers the range 10-316 eV. Notice 
that under the assumption adopted, the intrinsic 
number of photons emitted scales as VR, while 
the temperature as 1/ \fR. Therefore the highest 
intrinsic photon emissivity occurs at the lowest 
temperatures. The relation between hardness ra- 
tio and count rates for constant luminosity and 
varying area is shown in Fig. 9.c. The two curves 
corresponds to the same coronal parameters used 
in Fig. 9. a and 9.b. The same discrete points 
plotted in Fig. 9. a and 9.b are also marked, cor- 
responding to a decrease of R of 30% (and hence 
resulting in an increase of Tbb of 20%) with re- 
spect to the preceding one. 

The behaviour of HR is similar to that seen in 
the previous section. As long as UTbb ~ 60 eV 



the soft counts increase, while the hard counts 
are much less variable. This regime corresponds 
to the decreasing branches of the curves in Fig. 
9.c. Although HR is anticorrelated to the to- 
tal count rate independently of the coronal pa- 
rameters, the hard and soft count rates can be 
both correlated or anticorrelated, depending on 
whether the hard power law is steep or flat, re- 
spectively. 

The behaviour of HR vs. the total count rate 
is very different when VTbb ~ 60 eV. Since the 
total luminosity is constant, the total count rate 
is basically constant when the peak of the black 
body distribution is within the ROSAT band. 
This regime corresponds to the almost vertical 
branches of the curves in Fig. 9.c: HR dramat- 
ically increases at almost constant total count 
rate. In this regime the peak of the black body 
spectrum moves from the soft band to the hard 
band. The net result is that the two count rates 
are anticorrelated: as Tbb increases C^ 9 _ 2 ] in- 
creases while C[o.i_o.4] decreases. 

If the corona is pair dominated, variations of 
R should not happen with constant coronal pa- 
rameters. A reduction in the radius by a factor 
~ 10 (which roughly corresponds to four succes- 
sive points counting from left to right along the 
curves in Fig. 9.c) causes an increase in the opac- 
ity and hence a reduction in the coronal tempera- 
ture. For kTsB ~ 60 eV the variations of HR will 
be smaller with respect to cases with constant r 
and O. For larger values of VTbb the situation 
is similar to what seen before: a sharp increase 
of the hardness ratio with a constant count rate. 

In the case of a low pair density corona de- 
scribed in §3.2, the coronal parameters can vary 
without any change of size or luminosity of the 
soft component. In terms of ROSAT count rates, 
this is represented by the dashed lines in Fig. 
9.c. Those lines connect points relative to spectra 
with different coronal parameters, but same Tbb 
and luminosity. If for example an increase of r 
occurs, the ROSAT count rate increases as long 
as VTbb ~ 60 eV, slightly decreases for higher 
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black body temperatures (Fig. 9.c). The hard- 
ness ratio shows a slightly simpler behaviour: for 
low temperatures decreases, for kTsB ~ 30 eV it 
increases tending to become constant at higher 
temperatures. 

5. Spectral and Flux Correlations: Com- 
parison with Observations 

In the following we compare the results dis- 
cussed in the present sections with some of the 
current observations of soft-hard X-ray variabil- 
ity in Seyfert galaxies. 

The two recently studied cases of NGC5548 
and Mrk 766 observed with ROSAT and/or ASCA 
provide the most detailed available information 
on spectral and flux variability of Seyfert I galax- 
ies in the X-rays to date, and also seem to indi- 
cate opposite behaviours. We should warn that 
Mrk 766 is a rather extreme case, not represen- 
tative of the Seyfert I class. Indeed this source 
is classified Narrow Line Seyfert 1 (NLS1). 
Sources of this subclass tend to show rapid X-ray 
variability, unusual steep ROSAT spectra, and 
strong Fe II emission (Pounds & Brandt 1996). 

In the following we compare those observa- 
tions with the analysis developed in the previ- 
ous sections. We do not aim to "fit" any data, 
but simply to interpret (whenever possible) the 
general variability features in the context of the 
reprocessing scenario. 

5.1. NGC 5548 

The soft X-ray variability of NGC 5548 ob- 
served by ROSAT has been deeply analyzed by 
Done et al. (1995). The soft and hard counts 
seem to go up and down together, but not in 
a linear way. In general the source softens (in 
the ROSAT band) as it brightens. The hardness 
ratio roughly doubles (from 0.03 to 0.06) for a 
count reduction by a factor 3. Also a state where 
only the counts in the ROSAT soft band varied 
was observed. No simultaneous observations in 
the hard X-rays are available. Old GINGA data 



showed a power law index close to the canoni- 
cal value 0.9 plus a reflection hump (Nandra & 
Pounds 1994). 

The source requires a highly variable soft com- 
ponent, and a more steady hard one. This can 
be roughly interpreted assuming that variability 
in the soft band is due to variations in size of the 
emitting area, as in the case discussed in §4.1.2 
(see Fig. 9.c). The temperature of the ther- 
mal component should be less than 60-70 eV, 
not inconsistent with that found by Done et al. 
(1995). The intrinsic luminosity of the source 
should be roughly constant. We may notice that 
the lower value of the observed hardness ratio 
(0.03) is lower that the minimum value consis- 
tent with the reprocessing hypothesis (see Fig. 
9). The discrepancy may be due to the value of 
Nh for this source ( Nh = 1-65 x 10 20 cm -2 ; 
Done et al. (1995) ) which is lower than that 
used throughout our calculations (Nh = 2 x 10 20 
cm" 2 ), but it may also indicate that part of the 
reprocessed thermal flux does not contribute to 
the cooling of the corona, or alternatively, as sug- 
gested by Done et al. (1995), that there is an- 
other soft component unrelated to the hard one, 
possibly much broader than a single temperature 
black body. 

5.2. Mrk 766 

Mrk 766 has been observed by ROSAT and 
ASCA, providing a full coverage of the X-ray 
domain from 0.1 to 10 keV. The ROSAT data 
has been analyzed by Molendi, Maccacaro &: 
Schaeidt (1993), Molendi and Maccacaro (1994) 
and recently by Leighly et al. (1996) who also an- 
alyzed the ASCA data. ROSAT data show that 
the source hardens as it brightens, and can be 
roughly described by an almost steady soft com- 
ponent plus a variable hard one. In the ASCA 
band the source shows evidence of spectral vari- 
ability above 1 keV (Act ~ 0.4); the steepest 
state is brighter than the flat state in the ASCA 
band by a factor of ~ 2. The power law pivots at 
an energy close to 10 keV. If the power law ex- 
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tends up to few hundreds of keV, then the hard 
(and weaker in the ASCA band) state is actually 
more luminous than the soft state. 

Observations indicate that a[i_io] varies be- 
tween 0.6 and 1, while the total luminosity changes 
by factor of 2 at most. This indicates that either 
pairs are not dominant or, if pairs dominate, that 
the linear size of the emitting region is varying. 
In the latter case the required change in com- 
pactness is roughly a factor 10. Since the overall 
X-ray luminosity of the hard state is estimated 
to be 2 times larger than that of the soft state 
(see below), the region responsible for the soft 
spectrum should be 20 times smaller than that 
responsible for the hard one. The soft state, al- 
though fainter by a factor 2, should then be char- 
acterized by a Tbb roughly 3-4 times larger than 
that relative to the hard state. Data are consis- 
tent with a reduction of the excess normalization 
by a factor 2 as the power law steepens, but an 
increase of 3-4 times of Tbb in the soft state is 
not supported by the data. 

We therefore conclude that the spectral vari- 
ability in Mrk 766 could be due to changes of 
the scattering opacity of the corona without large 
variations of the intrinsic luminosity. The opac- 
ity of the corona is probably not dominated by 
electron-positron pairs. A variation of r from 
~ 0.3 - 0.4 to ~ 0.1 - 0.2 can explain the hard- 
ening of the power law, and possibly (part of) 
the decrease in the ASCA count rate. We pre- 
dict that the hard state is characterized by a 
coronal temperature ~ 150 — 200 kev (and hence 
an exponential cut-off with an e-folding energy 
Eq — 300—400 keV should be present), while the 
soft state should have kT ~ 80 — 100 keV (cor- 
responding to E c ~ 150 - 200 keV). The total 
X-ray luminosity of the hard state would then 
be 2 times larger than that of the soft state. 

Finally we note that the value of hardness ra- 
tio in the ROSAT band computed by Molendi 
& Maccacaro (1994) ranges from ~ 0.5 to ~ 1. 
According to Fig. 9. a this would imply a tem- 
perature of the soft component of ~ 130 eV, in 



agreement with the observed hardening with in- 
creasing intensity in the ROSAT band and with 
the ASCA results. 

6. Summary and Discussion 
6.1. Summary 

We have examined the expected spectral vari- 
ations within a model involving a hot corona 
emitting medium to hard X-rays by Comptoniza- 
tion, coupled to a cooler optically thick layer 
which i) provides the input of soft photons for 
Comptonization by the hot corona ii) intercepts 
a fixed fraction of the Comptonized photons and 
reprocesses them into soft thermal photons. The 
fraction of hard photons which are reprocessed 
and the fraction of soft photons which are Comp- 
tonized have been held fixed to the standard val- 
ues of 1/2 and 1 respectively which refer to a 
flattened, sandwich-like geometry. These values 
may be somewhat different in different sources 
but could be constant on "short" timescales, over 
which the structure of the emitting region is not 
expected to vary substantially. 

Even with these restrictions the spectral vari- 
ability is rather complex and depends on two 
main issues, the importance of pairs as a source 
of opacity and the expansion/ contraction of the 
active regions in the corona possibly associated 
with luminosity variations. 

The main point which emerges is that intrin- 
sic spectral variations are governed by changes 
in opacity (steeper spectra for larger r as long 
as t^I). For fixed L c /L s this implies that a 
steepening of a should correspond to a decrease 
of the temperature. This is a robust prediction 
of Compton cooled models with fixed geometry 
and could be tested by the recently launched 
XTE and SAX satellites. Should an opposite 
behaviour be found, the concept of a quasista- 
tionary corona above an accretion disk should be 
substantially revised in the sense of a much more 
chaotic situation. 

The relation between spectral shape and in- 
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tensity is univocal in the case of pair dominated 
coronae. In this case large changes in compact- 
ness (luminosity and /or size) are required for 
modest spectral variations. These should appear 
as either large changes in the 2-10 keV intensity 
or large changes in the temperature of the soft 
reprocessed component. The case of Mrk 766 
(which may be anomalous among Seyfert 1 galax- 
ies) seems to exclude a pair dominated corona. 

In the case of negligible pairs the relation be- 
tween luminosity and opacity is essentially un- 
known. Opacity variations could occur with lit- 
tle variations in luminosity and give rise to a 
pivoting of the medium energy spectrum which 
can yield little correlation and even anticorrela- 
tion between the hard and medium X-ray fluxes. 
For constant luminosity and flat spectral indices 
(«[2_io] ~ 1) the hard photon flux varies less 
and is anticorrelated with the photon flux in the 
medium energy band. 

The integrated luminosity in the soft repro- 
cessed emission is always proportional to the 
Comptonized luminosity. However its contri- 
bution to the counts in the band between the 
galactic absorption cut off and the Carbon edge 
depends strongly on the temperature which in- 
creases with increasing luminosity (for fixed size). 
This can give rise to anomalous observed be- 
haviours with larger variations in the soft band 
than in the medium one. 

6.2. Discussion 

From a more observational point of view, it 
may be useful to discuss in more detail what can 
be deduced about the intrinsic properties of the 
source from a given set of " observables" To this 
end we indicate the [0.1-2] keV band as the "soft 
S band", the [2-10] keV band as the "medium M 
band" and the band above 30 keV as the "hard 
H band". We define four observable quantities 
describing the continuum variability of a source 
in the X-ray domain: variation of the [2-10] keV 
spectral index Aa, variation of the [0.1-2] keV 
count rate AS, variation of the [2-10] keV count 



rate AM, and variation of the hard count rate 
AH. The four quantities defined are in principle 
independent, and hence all the possible variabil- 
ity behaviours can be described by a combination 
of possibilities. We will examine the most inter- 
esting ones. 

The main relevant property of an observed X- 
ray light curve is whether flux variability in the 
M band occurs keeping a constant spectral in- 
dex, or instead if it is accompanied to spectral 
variations. The two possibilities, together with 
simultaneous observations in the S band, lead to 
very different conclusions regarding the very na- 
ture of the observed variability. 

As discussed in §3 spectral variations in the 
medium band require a negligible amount of pairs 
in the corona, or, for large compactnesses, varia- 
tions of the source linear size R. Observations in 
the soft band can discriminate between the two 
possibilities, since the latter case would imply 
variations of Tbb and hence variations of the soft 
counts. Hence in case AS* ~ the intrinsic lumi- 
nosity of the X-ray source and the dimensions of 
the reprocessing region must not vary; the hard 
variability in flux and spectral index is due to 
variations of optical depth of the corona, with 
a negligible contribution of pairs. The model 
predicts that ap-io]; whenever steeper than 1, 
is anticorrelated to the medium count rate. A 
positive correlation holds if instead ap-io] ~ 1- 
The temperature of the corona decreases for in- 
creasing spectral index, and the flux in the H 
band is correlated to that in medium band if 
a [2-io] ~ 1) anticorrelated if ccp-io] ~ 1- The 
case with AS ^ is similar, but the observed 
variability properties can be due both to vari- 
ations of the intrinsic luminosity of the corona, 
and/or to variations of the linear size of the re- 
processing region. If the corona is pair domi- 
nated, the model predicts that the soft compo- 
nent temperature is correlated to ap-io] and an- 
ticorrelated to O. In case pairs are unimportant, 
it is not possible to draw any predictable relation 
between ap-io] and the medium or soft count 
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rates. A full coverage all through the X-rays to 
the 7— rays can determine in principle if the in- 
trinsic luminosity of the source varies. In case of 
negligible luminosity variations, changes in Tbb 
must occur to produce a non null AS. On the 
other hand, if intrinsic luminosity variations were 
detected, they alone could give raise to the ob- 
served variability in the soft band without asso- 
ciated variations oiTsB- In principle the nature 
of variability in the soft band (variations of in- 
trinsic flux or variations of temperature) can be 
tested by high spectral resolution observations. 

Different conclusions can be drawn if the vari- 
ability in the medium band occurs with a con- 
stant spectral index. The main result is that 
the intrinsic luminosity of the source must vary, 
but not the properties of the hot corona. This 
can happen in a pair dominated source if the 
variations of compactness are 2. An observa- 
tion of AM ^ 10 with a constant spectral index 
would imply that some mechanism different from 
pair balance works to keep the coronal optical 
depth constant. In any case, the model predicts 
that AS / 0. The total counts in the soft and 
hard bands are positively correlated, and much 
stronger variations occurs in the soft band than 
in the medium one. should not vary notice- 
ably, while the flux in H band follows the varia- 
tions of the M band. If A5 ~ the model can be 
ruled out, since any variation of the hard power 
law with a constant spectral index must be ac- 
companied by an observable variation in the lu- 
minosity of the soft component. If this does not 
occur, the soft component cannot be due to local 
reprocessing of part of the hard X-rays. 

It should be mentioned that, though the model 
considered in the paper constrains the ratio L c /L s 
to be constant ~ 2, in principle different be- 
haviours are possible if instead the ratio L c /L s is 
allowed to vary. This can occur if the height-to- 
radius ratio of the scattering region(s) is of order 
unity or larger (see, e.g. , HMG and Stern et al. 
1995). GH94 showed that in the case of a pair 
dominated homogeneous plasma with B ^ 1, an 



increase of L c /L s (and hence a decrease of a) 
would result in a lower temperature (and hence 
a relevant increase in r must occur) contrary to 
the case of r-driven spectral variations with con- 
stant L c /L s . Thus, at least in principle, broad 
band observations could be used to check the sta- 
bility of the Compton-to-soft ratio during spec- 
tral variations, and hence the role played by re- 
processing. 

Finally we wish to mention one aspect of the 
model neglected so far, namely time lags be- 
tween different bands. In the most basic picture 
of Comptonization models, the hard X-rays re- 
spond to variations of the soft photon flux, and 
hence variations in the hard band should fol- 
low variations of the soft thermal photons, ex- 
hibiting also a steeper power spectrum (see e.g. 
the discussion in Papadakis & Lawrence 1995). 
Data in this field are sparse and often contradic- 
tory (see Tagliaferri et al. 1996 and references 
therein). EXOSAT data of NGC 4051 seem in- 
dicating that more power at high frequency is 
radiated by means of high energy photons (Pa- 
padakis & Lawrence 1995), contrary to what ex- 
pected from the basic Comptonization theory. 

The above considerations assume that spec- 
tral variability is due to variations of the primary 
soft photon input. In the model under investiga- 
tion here however the soft photons are supposed 
to be generated by mere reprocessing of a frac- 
tion (~ 50%) of the primary hard X-rays. Varia- 
tions of the physical conditions occur primarily in 
the corona. For example, we may expect that a 
rapid variation of the heating rate appears in the 
Comptonized photons roughly at the same time 
at every energy, and eventually in the soft com- 
ponent. The situation appears to be extremely 
complicated, and we believe that data interpre- 
tation can not avoid to deal with such complica- 
tions (see, e.g., Nowak &: Vaughan 1996). In view 
of forthcoming better timing data, in a work in 
progress we aim to perform a detailed analysis of 
the timing properties of an active Comptonizing 
corona. 
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Fig. 1. — The total (direct power law + Comp- 
tonized reflection) spectral index in the [2-10] 
keV range ap-io] is plotted against the coronal 
scattering optical depth r. The open symbols 
connected with dashed lines refer to spectra com- 
puted with kTsB = 100 eV. The closed symbols 
connected with solid lines are for JzTbb = 50 eV. 
The squares indicate that a face-on line of sight 
is assumed, the triangles a 60° inclination. 
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Fig. 2. — The total (direct power law + Comp- 
tonized reflection) spectral index in the [2-10] 
keV range ap-io] is plotted against the coronal 
electron temperature kT. Symbols as in Fig. 1. 
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Fig. 3. — The difference between the spectral in- 
dex of the [2-10] keV intrinsic continuum and 
that of the [2-10] keV total spectrum (contin- 
uum +reflection) is plotted vs. the optical depth 
r. Symbols as in Fig. 1. 



Fig. 4. — The total spectral index ap-io] is plot- 
ted vs. the photon flux in the [2-10] keV band 
(arbitrary normalized) in the case of a pair dom- 
inated corona. Squares connected with solid line 
refer to a case with a face-on line of sight, trian- 
gles connected by dashed line to an inclination of 
60°. kTsB = 100 eV in both cases. 
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Fig. 5.— The photon flux in the [30-100] keV 
band is plotted vs. the photon flux in the [2- 
10] keV band in the case of a pair dominated 
corona. Normalization is arbitrary. The curves 
are for kTsB = 100 eV, and squares indicate a 
face-on line of sight, triangles a 60° inclination. 
The value of the total spectral index ap-io] m 
the [2-10] keV range is indicated for some of the 
spectra. 



Fig. 6. — Comptonized spectra for different 
values of (r, 0) equilibrium values [(0.63,0.11), 
(0.32,0.21), (0.2,0.3), (0.1,0.5)]. kT BB is set to 
100 eV, and a face-on line of sight is considered. 
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Fig. 7. — Same as Fig. 4, but in the case of 
r variations with constant luminosity. Squares 
with solid line are for a face-on line of sight, tri- 
angles with dashed line for a 60° inclination. 
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Fig. 8. — Same as Fig. 5, but in the case of 
r variations with constant luminosity. Squares 
with solid line indicate a face-on line of sight, 
triangles with dashed line a 60° inclination. 
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Fig. 9. — (a) The hardness ratio in the ROSAT 
band plotted vs. the temperature of the repro- 
cessed component. Squares refer to spectra com- 
puted with r = 0.1, = 0.5, while triangles 
with r = 0.63, = 0.11. (b) The hardness ratio 
is plotted vs. the total count rate in the ROSAT 
band for a case of pure luminosity variations. 
Squares and triangles same as (a). Dashed lines 
connect spectra having equal kTsB- ( c ) Same as 
(b) but for pure variations of the emitting area. 
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